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We show that it is possible to achieve maximally entangled mixed states of two qubits from the 
singlet state via the action of local non-trace-preserving quantum channels. Moreover, we present a 
simple, feasible linear optical implementation of one of such channels. 
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I. INTRODUCTION 

In a recent paper, Ziman and Buzek have demon- 
strated that it is impossible to transform the singlet 
state of two qubits in a maximally entangled mixed state 
(MEMS), via local maps £®I m. Such maps describe 
the action of quantum channels Cs acting on a single 
qubit of the initial singlet state. When a channel is local, 
that is when it acts on a single qubit, the corresponding 
map is subjected to some restrictions. This can be eas- 
ily understood in the following way: Let Alice and Bob 
be two spatially separated observer who can make mea- 
surements on qubits a and b, respectively, and let pin and 
Pout denote the density matrices describing the two-qubit 
quantum state before and after the channel, respectively. 
In absence of any causal connection between Alice and 
Bob, special relativity demands that Bob cannot detect 
via any type of local measurement the presence of the 
channel Cs in the path of photon a. Since the physics of 
qubit b is described to Bob by the reduced density matrix 
TrpoutU, the locality constraint can be written as 



Pout 



Pout ~ Pu 



If we write explicitly the map £ (g) X as 

Pin ^ Pout = ^ (A^ /) Pin {Al (g) /) , 

then Eq. |^ becomes 

'^iPin)H,k3(^'^AlAf,'j = ^(Pi„)fcj,fcj, 
k.l a k 



(1) 



(2) 



(3) 



which implies the trace-preserving condition on the local 
map £ (g) X: 



(4) 



Local maps that do not satisfy Eq. Q are classified as 
non-physical, and are not investigated in Ref. [J. 

In this paper we show that under certain circum- 
stances, it may be meaningful to consider the action of 
non-trace-preserving maps, as well. In particular, we give 
two simple examples of local non-trace-preserving maps 



that generate maximally entangled mixed states of two 
qubits from the singlet state. Two-qubit MEMS states 
may exist in two subclasses usually denoted as MEMS 
I and MEMS II In Section II, we furnish an ex- 
plicit representation for two maps M. and /C that gen- 
erate MEMS I and II states, respectively. In Sec. Ill 
a feasible linear optical implementation of the quantum 
channel Cm corresponding to the map M. is given. In 
Sec. IV we introduce an all-unitary linear optical model 
for Cm and, via a rigorous QED treatment, we show how 
the "non-physical" map AA arises in a natural manner. 
Finally, we draw our conclusions in Sec. V. 



II. NON-TRACE-PRESERVING MAPS 

In this section we introduce two non-trace-preserving 
maps M and JC that generate MEMS I and II states, 
respectively, from an initial singlet state of two qubits. 



A. MEMS I map M 

Let Pin represent the initial state of a single qubit that 
is transformed under the action of the map M. as: pin ^— *■ 
Pout, where 



Pout = E M^pinMj; , 



(5) 



and Mo = = Ml, 



M2 = V2(l-P) 



1 
0/' 



Ms = VP K 



-1 




(6) 

where 2/3 < p < 1. This map is not trace-preserving nor 
unital, since 

^MtM, = ^M,Mt = (2^P 0^1, (7) 



M=0 



^1=0 



where I„ denotes the n x n identity matrix. The apparent 
non-physical nature of this map can be displayed if we 
use the Pauli matrices to rewrite 



Mo 



l2+cr^), mi = -i^(Ty, (8) 



2 



and substitute Eq. © into Eq. ^ to obtain 



where Ki = 0, and 



Pout = 2 [(-"^ ~ P)'^'" ^ <^yPin(Ty + (1 - p) CTzPinfTz 
+ (1 --P){An,0-^}], 

(9) 

where the anti-commutator term ({a, 6} = a6 + 6a) is 
clearly responsible for non conservation of the trace. 

Now, let us consider the map Ai as representative of 
the local quantum channel Cm acting on a single qubit 
belonging to an entangled pair prepared in the initial 
state Pin (note that now pin denotes a two-qubit state, 
therefore it is represented by a 4 x 4 matrix). The two- 
qubit map M can be written as 



Kn = X - 



1 




1 

1 



-1 

1 



(15a) 



(15b) 



(15c) 



and < p < 2/3. As in the case of A4, this map is not 
trace-preserving nor unital, since 



Pout = ^N^PinNj,, 

where = (g) I2, namely No = = Ni, and 



(10) 



No 



10 




(lla) 



N3 = VP 



/O -1 
-1 
10 

\0 1 



(lib) 



Let 10") = (|01)-|10))/V2 be the two-qubit input singlet 
state represented by the density matrix ps'. 



/O 0\ 

i -i 

Ps=r)(^-i= 4 r^ ° 

\ 0/ 

A straightforward calculation shows that 



(12) 



PI 



1-1=0 



f f 

1-p 


£ 

2 



(13) 



f 



which represent a MEMS I state. 



B. MEMS II map IC 

As before, let pin represent the initial state of a single 
qubit that transforms under the action of the map K. as: 
Pin 1-^ Pout, where 



Pout 



(14) 



M=0 



fj.=0 



2 
1 



(16) 



A straightforward calculation shows that the two-qubit 
map IC realized by = (g) I2 produces MEMS II 
states when acting upon the singlet state H12|l: 



Pii 



3 

fj.=0 



i\ 












f i 



(17) 



where Li = 0, and 



-0 = -\ 77 



10 

10 







(18a) 



L2 = A / - - - 



10 

P I 1 

2 I 1 

10 



(18b) 



/O -1 
-1 
10 

Vo 1 



(18c) 



III. 



LINEAR OPTICAL IMPLEMENTATION OF 
THE CHANNEL Cm 



The layout of the experiment we propose to create 
MEMS I states is illustrated schematically in Fig. ^ 
Two photons in the singlet state {\HV) - \VH))/V2 
emerge from the down-converter. Here H and V are 
labels for horizontally and vertically polarized photons, 
respectively. Photon b goes directly to detector Df,, while 
photon a goes to C^t and then to detector D^. Cj^ is a 
linear optical two-port device that is illustrated in detail 



3 



matrix B as: 



photon a 



Cm 



laser 
pump" 



down 
converter 



photon b 



FIG. 1: Sketch of the proposed experimental setup. The 
box Cm represents the quantum channel. A photon from an 
intense laser pump is split into the pair (a, b) by the down- 
converter. Detectors Da and Dj, permit a tomographically 
complete reconstruction of the two-photon quantum state. 
Further details are given in the text. 



in Fig. 121 Supposedly, detector Da does not distinguish 
which output port of Cm the photon comes from: This 
is our mechanism to induce decoherence. Photon a en- 
ters port 1 and can be either transmitted to path 1 or 
reflected to path 2 by the 50/50 beam splitter BS; vac- 
uum enters port 2. Let the square bracket vectors ' 



and 







represent the two orthogonal position states 



(or, spatial modes) of a photon travelling in paths 1 and 
2, respectively. Analogously, let the parenthesis vectors 

Q ^ and {^\^ represent the two polarization states of 

a photon polarized along a horizontal and a vertical di- 
rection, respectively. Described in these terms, the 50/50 
beam splitter performs a linear transformation restricted 
to the mode space only; it can be represented by the 2x2 



photon a 
(port 1) 



A Vh 



/ 



BS 




vacuum 
(port 2) 

FIG. 2: Detailed scheme of the quantum channel Cm- BS 
denotes a 50/50 beam splitter, A is a beam attenuator, Ph is 
a linear polarizer that selects horizontally polarized photons, 
and R is a polarization rotator oriented at = 7r/2. 



1 i 
i 1 



(19) 



where the relative phase shift of 7r/2 between the 
transmitted and reflected amplitudes ensures unitarity: 
BBt = I2. The attenuator A can be simply represented 
by a scalar function exp(— a), where a > 0. The linear 
polarizer P// performs a linear transformation restricted 
to the polarization space only. It can be represented by 
the projection matrix H as: 



H 



1 




(20) 



Finally, the polarization rotator R can be represented by 
the orthogonal matrix R in the polarization space as: 



R 



cos U — sm f 
sin 6 cos 6 



-1 

1 



(21) 



If with P and Q we denote the two complementary mode- 
space projectors 



(22) 



then the total 4x4 transmission matrix T representing 
Cm can be written as: 



"10' 


Q = 


■ 


0" 








1 



T = R(»P.B-f e""H® Q.B, 



(23) 



where the low dot " . " denotes the ordinary matrix prod- 
uct. 

Let I in) be the quantum state of a photon entering Cm 
through port 1: 



(24) 





' 1 " 








1. In these terms, the two-mode 



where 10^^ -f \(t>v\ 
output state |out) leaving Cm can be written as 

lout) = Tlin) 



R|(^ 

1 
71 



.P.B|?/') 



-0y 

4>H 



Q.B| 



(pH 







1 

(25) 



An elementary calculation shows that 
I out) (out I ^ T I in) (in I 

= R|(/))(0|R^ (g)P.B|i/;)(?/>|B^p1' 

+e-2"H|0)((/)|H^ ® Q.B|?/')(V'|B1'.Q^ 

-H{e-"R|0)((/)|H^ ® P.B|7/;)(V'|Bt.Q^ 

+H.C.}, 



(26) 
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where H.c. stands for Hermitian conjugate. 

Since, by hypothesis, detector does not distinguish 
a photon exiting port 1 from a photon exiting port 2, the 
detected output state can be obtained from | out) (out | by 
tracing over the detected but unresolved position states: 



Tr[|out)(out|] = i (R|0)(^|R^ + e-2"H|0)(0|Ht) , 

(27) 

where triviaUy Tr[P.B|7/>)(V'|Bt.Qt] = 0, and 
Tr[P.B|7/')(V'|B^Pt] = 1/2 = Tr[Q.B|?/')(i/'|Bt.Qt]. If 

we define Jout = Tr|" |out)(out|] and pin = \(j)){4>\, then we 
can rewrite Eq. (|27|l as 



i / + ve 



MgPinM^ ) , 



(28) 



where Eqs. 10, H2U|) . and H21|l have been used. If we 
choose a — a{p) such that 



Pe"'" , ^ . ^ 1, 2(1 -p) 

— - = 1 alp) = m 

2(1 -p) 2 



then Eq. H28|l can be rewritten as 



1 ^ 1 



M=0 



2p 



P 



Pout, 



(29) 



(30) 



where Eq. (0 has been used. Note that a{p) > for 
2/3 < p < 1, as expected for an attenuator. Equation 
()30|1 shows that the scheme shown in Fig. |21 actually 
implements the map Ai. Moreover, from Eqs. H24l28t 
ESJl it follows that 



Tr (Jout) 



2(1 -p) 



\<I>h\ 



(31) 



< Tr(Jout) < 1, 



for 2/3 < p < 1 and < \4ih\ < 1. This means that even 
in the worst case (p = 1) there is still a 50% of probability 
to detect a photon in our scheme. 



IV. RIGOROUS QED TREATMENT 

It was pointed out |3| that the map Ai corresponds 
to a non-physical quantum channel Cm- Conversely, in 
the previous section we have shown that a physical linear 
optical implementation of Cm is actually feasible. The 
resolution of this apparent paradox lies in the concep- 
tual difference that exists between the "quantum state 
of two qubits" , and the "measured quantum state of two 
qubits" . The latter can be reconstructed only after Alice 



and Bob have performed coincidence measurements 
that is only after they have established a communica- 
tion and have compared their own experimental results. 
Therefore, a measured MEMS state generated by a local 
channel does not raise any causality issue. In this spirit, 
we will soon show how the map A4 can be derived from 
an all-unitary model for the channel Cm (Fig- 3)- Such a 
unitary channel reduces to a non-unitary one when Alice 
restricts her measurements to two output ports only (1 
and 2), leaving the other two (3 and 4) undetected. How- 
ever, note that in principle Alice could use an additional 
detector D^^^' coupled to ports 3 and 4 to generate a 
"conditional" MEMS state: When a photon pair is cre- 
ated by the down-converter and detector Di^**^ does not 
fire, then a conditional MEMS state is being transmitted 
through the channel. 

Let indicate with and 6q the annihilation operators 
of photons a and 6, respectively. Greek indexes a, (3, . . . S 
{0, 1} label polarization modes of the field, while Latin 
indexes i,j,... e {1,2,3,4} label spatial modes of the 
field. The latter modes represent the four paths shown 
in Fig. 3. Described in these terms, the two-photon input 



3 

PBS * 



4 

' A 



photon a 
(port 1) 



BS 




R 



vacuum 
(port 2) 

FIG. 3: All-unitary linear optical realization of the quantum 
channel Cm ■ PBS is a polarizing beam splitter that transmits 
only horizontally polarized photons. The beam attenuator A 
is implemented by a variable-reflectivity beam splitter. The 
two additional modes 3 and 4 ensure the whole unitary nature 
of Cm- 

singlet state can be written as 

\m) = ^[a\„b\-albl)\0), (32) 

where |0) denotes the vacuum state. Each linear optical 
element present in the quantum channel shown in Fig. 3, 
can be represented by a unitary operator U that evolves 
the state vector \tb) as 



and the operator X either as 

X ^ U^XU, 



or as 



X >-» uxu'f. 



(33) 
(34) 
(35) 
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In particular, if the annihilation operator evolves as 

4 1 

ttia U^'UiaU = ^^^^Siajf3ajf3, (36) 

then it is easy to see that 

4 1 

0-L Ual^U'' ^J2Y1 '^i/3.»""]/3- 



(37) 



Within this formalism, the beam splitter BS is described 
by the matrix 

Sjl3,ia — ^jiSpa (38) 

where B is explicitly given in Eq. (|19|l . This leads to the 
field operators transformation 



that modifies the input state |in) to: 
The effect of the rotator R is very simple: 



(39) 

(40) 
(41) 
(42) 



Next, the polarizing beam splitter PBS can be described 
by a 4 X 4 unitary matrix that couples both spatial and 
polarization modes, as 



"10 " "11' "11 " ^"10' 

and it changes the two-photon states to 



f 4o\ / 1 \ / aL \ 



21 

''30 



V4i/ 



i 
10 
Vo i 



21 

-'so 



(43) 



V4i/ 



As a result of this transformation, the two photon state 
after the PBS can be written as 

|in) ^ J (alAi + ^4A + alAo + 4Ao) |0). (44) 



Finally, the attenuator A can be described in a unitary 
fashion by modelling it as a variable-refiectivity beam 
splitter such that: 



^2a 



(45) 



where < T < 1, and + i?^ ^ 1. This last optical 
element produces the output state |out), where 



but 



(a\Ao-Ra\Ai)] |0) 



+ 

IV'12) + IV'34). 



(46) 



In Eq. denotes the two-photon state restricted 

to the pair of modes and ('0i2|'034) = 0- Since each 

transformation performed by each linear optical element 
present in the quantum channel is unitary, the output 
state |out) is still normalized: (outjout) = 1. Now we 
can trace over the spatial degrees of freedom in the usual 
way obtaining: 



p = Tr [|out)(out|] = P12 + P34, 



where 



P12 



/ 1 1 

T2 



V 1 1 



(47) 



(48) 



and 



P34 



1 

4 















1 
















1 


















(49) 



The total density matrix p has trace equal to one, while 
the truncated density matrices pi2 and P34 have nonunit 
trace: 



Tr(pi2) = i(l 



1 

2^' 



(50a) 



Tr(p34) = ^(l-TV2)^l-^. (50b) 

This simple result shows that each truncated density ma- 
trix cannot be generated by a trace-preserving map. In 
particular, we easily recover our result Eq. H13|l by divid- 
ing Eq. by Eq. (|50a|l . This is the goal of the present 
section. 

To conclude, it may be instructive to calculate sepa- 
rately the reduced density matrices Trpyj^, obtained by 
tracing over the degrees of freedom of photon /, where 
(i,j)e{(l,2),(3,4)}and/ = a,6: 



1/1 



Trpi2|„ =4(0 1 + ) ' '^'■^121;, 



1/1 



034 If, 



1 / 1 + r2 
1 

1 / 1 - r2 
1 

(51) 



From these results we learn that 



Trpi^L + Trp34L = 1(1 I ] = T^ftL , (52a) 



Trpi2|, + Trp34|, = ^ ( J J ) = Trps|, , (52b) 

that is, when all spatial modes of the two photons are 
properly accounted for, locality requirements are fully 
satisfied. 
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V. CONCLUSIONS 

Equations (O, (CZIl, CHll, and are the main 
results of our preliminary work on generation and mea- 
surement of maximally entangled mixed states. In this 
paper we have shown how it is possible to generate both 
MEMS I and II two-qubit states from the singlet state by 
using only local, non-trace-preserving quantum channels. 
Moreover, we provided for the scheme of a simple linear 
optical experimental setup for the generation of photonic 



MEMS I states. Such a scheme, which exploit spatial 
degrees of freedom of the photons to induce decoherence, 
is currently being tested in our laboratory. 
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